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_ (57) Abstract: The invention relates to a method for cancelling feedback in an acoustic system comprising a microphone, a signal 
^ path, a speaker and means for detecting presence of feedback between the speaker and the microphone, the method comprising 
q providing an LMS algorithm for processing the signal; where the LMS algorithm operates with a predetermined adaptation speed 
when feedback is not present; where the LMS algorithm operates an adaptation speed faster than the predetermined adaptation speed 
Q when feedback is present, and where the means for detecting the presence of feedback is used to control the adaptation speed selection 
^ of the LMS algorithm, where the feedback detection means comprises bandwidth detection means for determining the presence of a 
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TITLE 

Feedback cancellation using bandwidth detection 

5 TECHNICAL FIELD 

The present invention concerns hearing aids. In many hearing aids, for example 
"In-the-ear" (ITE) and "Behind-the-ear" (BTE), the microphone and the receiver 
(telephone) components are placed close to each other. This may result in that 
10 the sound produced by the receiver leaks back into the microphone. This may 

occur when the hearing aid shell or the ear mould does not fit sufficiently tight in 
the ear canal. Given enough amplification in the hearing aid, the loop gain of the 
system may exceed 0 dB at some frequency and a feedback oscillation may be 
produced. 



15 



BACKGROUND OF THE INVENTION 



The present invention is based on algorithms previously proposed in the 
literature. The invention concerns a number of algorithm modifications, which 
20 overcome some of the limitations of other systems used for feedback reduction in 
hearing aids. 

The invention relates to a feedback cancellation algorithm, which does not need 
an artificial noise signal in order to estimate the feedback transfer function. The 

25 input signal received from the environment, or the feedback oscillation signal, is 
used to drive the estimation process. In this fashion, the hearing aid user does not 
listen to an added noise signal, and a higher sound quality is possible. However, 
it is well known that such 'no-noise' algorithms can have audible side effects 
under certain circumstances, especially when environmental signals with long 

30 autocorrelation functions are present at the microphone. 



The autocorrelation function for a signal describes the average correlation 
between two signal values, which are separated by a time difference "Lag". In 
loose terms, the autocorrelation function describes how "predictable" a signal 
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value is, given the other samples in the signal. Some signals, for example periodic 
signals, are highly predictable and, correspondingly, the autocorrelation function 
does not vanish even for large values of Lag. Other signals, such as white noise, 
are very little predictable, and their autocorrelation function quickly vanish for 
5 increasing values of Lag. For signals with a long autocorrelation function, a future 
sample value can be predicted with a high degree of confidence, given the past 
samples. In other words, new samples of the signal do not provide much new 
information. Careful analysis of feedback cancellation systems reveal that signals 
with long autocorrelation may drive the adaptive system to produce poor 
10 estimates of the feedback path. 

It is the objective of the present invention to provide a method and a hearing aid 
for feedback cancellation, which improves the result of the feedback cancellation 
by having fewer audible side effects and thereby gives an improved user comfort. 

15 

SUMMARY OF THE INVENTION 

According to the invention the objective is achieved by a method, which 
comprises the feature of claim 1. According to the invention the first objective is 
20 likewise achieved by a hearing aid, which comprises the features of claim 6. 

Hereby a low adaptation speed, which generally improves the sound quality for 
signals with long autocorrelation functions, is applied when no feedback 
oscillation is present and a high adaptation speed, which is desirable to reduce 
25 feedback oscillations quickly, is applied when feedback oscillation is present 

hereby maintaining the preferred mode when feedback is not present and quickly 
changing the mode essentially without audible oscillations. This results in fewer 
audible side effects and an improved user comfort. 

30 Advantageous embodiments are described in the dependent claims 2-5 and 7. The 
function and effect of these is explained in connection with the preferred 
embodiment. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 is a schematic diagram showing a previously known system used for 
feedback cancellation; 

5 

FIG. 2 is a schematic diagram showing an embodiment of the system according to 
the present invention. 

FIG. 3 is a schematic diagram showing the feedback detection system according 
10 to the invention; 

DETAILED DESCRIPTION OF THE DRAWINGS 

A well-known principle for feedback cancellation in hearing aids is shown in fig. 
15 1. All the components described below, except blocks (1), (5) and (50), operate in 
the discrete time domain. 

The components are as follows: (1) is a microphone which picks up the sound 
from the environment (51) ("External input") and the feedback signal (52) 

20 ("FBSignal"); (2) is a microphone amplifier and an analog-to-digital converter 
(A/D); (3) is the hearing aid amplifier with filters, compressors, etc.; (4) is a 
digital-to-analog converter and a power amplifier; (5) is the hearing aid receiver; 
(50) is the acoustic feedback path (outside the hearing aid); (6) is a delay unit 
whose delay matches the delay through the components (4), (5), (50), (1) and (2). 

25 (7) is an N-tap finite impulse response (FIR) filter which is intended to simulate 
the combined impulse response of components (4), (5), (1), (2) and (50). (8) is an 
adaptive algorithm which will adjust the coefficients (9) of the filter (7) so as to 
minimize the power of the error signal (10). 

30 The algorithm (8) is well known as the Least Mean Square (LMS) algorithm. The 
algorithm requires a reference signal (11), which is used to excite the path 
consisting of the components (4), (5), (1), (2) and (50). The correlation between 
the reference signal (11) and the error signal (10) is used to compute the 
adjustment of the coefficients (9). 
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The system utilizes the output signal (11) from the hearing aid amplifier block (3) 
as a driving signal for the LMS algorithm, thereby eliminating the need for a 
disturbing noise in the receiver (5). 

5 

For some external input signals, the LMS based algorithm used in the application 
shown in fig. 1 is known to have difficulty adjusting the coefficients (9) as 
desired, i.e. to match the path consisting of components (4), (5), (1), (2) and (50). 
The difficulties are greatest for signals with long autocorrelation functions. 

10 Mismatched coefficients may lead to audible side effects, which can be very 
disturbing to a hearing aid user. Those may comprise audible oscillations and 
change in gain characteristics and frequency characteristics. One general remedy 
against this problem is to use a low adaptation speed, but this leads to poorer 
performance of the system because the coefficients cannot track changes in the 

15 acoustic feedback path (50) quickly, resulting in a long feedback cancellation 
time. 

The basic system shown in fig. 1 may be improved in various ways to minimize 
the side effects associated with certain input signals. Many authors have 
20 proposed additional system blocks, which will improve the sound quality while 
maintaining an acceptable adaptation speed. 

The present invention is based on the system diagram shown in fig. 1, and the 
invention consists of additional features, which will improve the sound quality 
25 and maintain an acceptable adaptation speed. 

FIG. 2 shows the block diagram of the general system and the components of the 
invention. 

30 The embodiment shown includes three features: Adaptation rate control, a 
frequency-selective adaptation procedure, and a feedback oscillation detector. 

Adaptation rate control 
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Two well known operation modes for the LMS algorithm are the "standard" mode 
and the "normalized" mode. In the "standard" mode, the coefficients are updated 
by an amount that depends on the short-term power of the error signal and the 
reference signal. This means that the update rate is faster when more powerful 
5 signals are processed by the hearing aid. In the "normalized" mode, the update 
rate is made nearly independent of the signal power, due to a normalization of 
the update equation. 

As described earlier, a low adaptation speed generally improves the sound 
10 quality for signals with long autocorrelation functions. In contrast, a high 
adaptation speed is desirable to reduce feedback oscillations quickly. 

Other authors have previously proposed changing the adaptation rate factor (often 
known as "ji") when feedback oscillations are detected. Although this does 
15 increase the adaptation speed, it also allows coefficients to deteriorate 

proportionally faster, in those situations where signals with long autocorrelation 
functions are present at the hearing aid input. 

In the present invention, the fact that feedback oscillations often have a high 
20 power is used. In many hearing aids, the output level is limited by compressor 

circuits, and in many cases the maximum output level is well above the normally 
used output level, for example when speech and other environmental signal are 
present. It is therefore assume that the feedback oscillations have a higher power 
than the environmental signal, in most cases where feedback problems exist. 

25 

Additionally, the feedback oscillation has the desirable property that its 
frequency is generally equal to the frequency where the loop gain currently is 
highest, i.e. where the fastest adaptation is needed. 

30 For the reasons mentioned above, it is very effective to utilize the feedback 
oscillation signal itself as a driving signal for the adaptation. 



When the "normalized" adaptation approach is used, the high-power feature of 
the feedback oscillation is not utilized. If, instead, the "standard" update 
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approach were used, the high power feature of the feedback oscillation would be 
utilized. At the same time, however, stronger signals in general would cause a 
higher adaptation speed, which could lead to more autocorrelation problems. 

The present invention introduces a new normalization scheme, which will 
generally maintain the low adaptation speed and the normalized operation mode, 
except when a feedback oscillation is detected. When a feedback oscillation is 
detected, the system is switched from normalized operation to standard operation 
by the switch (13), and the full power of the feedback oscillation signal is 
therefore allowed to adapt the coefficients. During "standard" operation, the 
update parameter (14) is chosen to such a value (53) that the external input (51) 
produces approximately the same update rate as it would in "normalized" 
operation. Assuming that the external input signal (51) maintains nearly constant 
properties before and during the feedback oscillation, the switch of normalization 
procedure will be nearly transparent to the external signal (51). This ensures that 
the sound quality remains high, even though the adaptation speed has been 
increased due to the higher power in the feedback oscillation. The update 
parameter (53) to be used during standard mode is estimated in component (12) 
before the feedback oscillation is detected. During intervals of feedback 
oscillations, controls signal (15) prevents (12) from updating the parameter (53). 

The switch from normalized mode to standard mode may be controlled by a 
feedback oscillation detector (49) through its output signal (15). The switch (13) 
may also be controlled by other conditions, which could result in feedback 
oscillations, for example when the acoustic feedback is rapidly decreased. 

The adaptive LMS algorithm (8) may be implemented as the following set of 
equations: 

Normalized operation: 



hk(n + 1) = hk(ri) + 



a-r(n-k)e(n) 



9 p = L.N 



p 



k=l..N 



(El) 
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hk(n + \) = hk(n) + 



ar(n-k)-e(n) 

b + LTsum 



(E2) 



5 In these equations, h k (n) is the k'th coefficient in the FIR filter at sample time n; a 
is a constant which determines the general adaptation speed of the algorithm 
(this constant is sometimes referred to as b is a small constant which 
prevents division by 0 for very small values of the reference signal; N is the 
number of coefficients in the filter (7); r(n) is the reference signal (30) sample 
10 value at time n; e(n) is the error signal (28) sample value at time n; and LT Sum is a 
value computed as described below. 

The sum term of the denominator of El is equal to the signal (54). LT 8um is equal to 
the signal (53), 



LT iUB1 (equal to (53)), which is computed by component (12), may be updated 
according to eq. (E3a): 



ot LT and p LT are time constants, which control the length of the exponential 
window over which the value of LT 8Uin is computed. 

Eq. (E3a) should not be updated while a feedback oscillation is present, since 
LT W1II should reflect the long-term value of SumSq for segments without 
oscillation. Once the feedback oscillation has disappeared, eq. (E3a) may be 
updated again. 



15 



LT*um(n + 1) = LTsum(n) • f5u + SumSq(n) • out 



(E3a) 



20 In equation (E3a) SumSq(n) is defined as follows (E3b): 




(E3b) 
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In El and E3b, the reference signal r(n) is used for normalizing the update 
equation. However, other signals in the system shown in fig. 2 may also be used 
instead of r(n). In the literature, the error signal e(n) has been used instead of r(n) 
5 for normalization; and even combinations of r(n) and e(n) have been used. The 
present invention will work for any type of normalization, in which the 
denominator in El and E2 is increased when the power level in the feedback loop 
consisting of (1), (2), (3), (4), (5) and (50) is increased. 

10 Frequency-selective adaptation 

Many feedback cancellation systems proposed earlier contain some form of 
frequency weighting of the signals which enter the LMS algorithm (8). The 
purpose of such weighting is to attenuate frequency ranges in which the 
15 autocorrelation of the external input signal (51) is long, and thereby reduce the 
possibility of poorly adjusted coefficients and poor sound quality. Several 
possibilities exist for frequency weighting. Various combinations of fixed and 
adaptive filters have been suggested in the past. 

20 In the present invention, steep highpass filters with high attenuation (20) are 
included in the inputs to the LMS algorithm. The purpose of these filters is to 
prevent low frequency contents from the reference signal (11) from entering the 
LMS algorithm. The cutoff frequency for the highpass filters (20) must be lower 
than the lowest frequency for which feedback cancellation should take place, and 

25 otherwise as high as possible. 

With the highpass filters (20) in place, the LMS algorithm (8) would not 
experience an increased level of the error signal (10) when the coefficients (9) are 
poorly adjusted in the low frequency range. Filter (7) with poorly adjusted 
30 coefficients, combined with components (3) and (6), may lead to a system with a 
high loop gain, and instabilities may result. 

In order to avoid this problem, a parallel feedback cancellation filter (21) is 
added. This filter is intended to provide low frequency information to the LMS 
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algorithm. The two filters (7) and (21) use identical coefficients (9). While filter 
(7) is designed to simulate the path consisting of components (4), (5), (1), (2) and 
(50), filter (21) is designed to simulate the artificial path (25) with an impulse 
response of constant '0'. The adder (33) computes an error signal as the difference 
5 between the desired '0' output and the actual output (34) from the filter (21). The 
error output (10) from the high frequency range and the error output (27) from the 
low frequency range are combined into a single error signal (28) which is fed to 
the error input of the LMS algorithm (8). In order to generate a low frequency 
signal as input to the filter (21) and to the reference input to the LMS algorithm, a 

10 noise generator (22) is included. The noise generator output (29) is lowpass 
filtered by a fixed filter (23). The cutoff frequency for the lowpass filter (23) is 
selected approximately equal to the cutoff frequency of the highpass filters (20), 
to obtain a reasonably flat input spectrum to the LMS algorithm. The low 
frequency signal (32) and the high frequency signal (31) are combined by the 

15 adder (24) to form the complete reference signal (30) for the LMS algorithm. 
Clearly, the components (25) and (33) may be removed immediately, and the 
signal (34) can be connected to the signal (27). 

The noise generator (22) may be implemented by randomly swapping the 
20 numerical sign of each sample of the signal (35). In other words, for each sample 
instant it is randomly decided whether the sample value should be multiplied by 
1 or by (-1). The advantage of using this type of noise generator is that noise 
samples at (35) and at (29) always have the same amplitude. The power spectrum 
of the reference signal (30) is therefore reasonably balanced at all times. In the 
25 literature, the noise generated as described above is sometimes referred to as 
'Schroeder' noise. 

Feedback oscillation detector 

30 Feedback oscillations may be produced by a system which contains an amplifier 
and a feedback loop, under some circumstances, A hearing aid with acoustic 
amplification, combined with an acoustic path from the hearing aid telephone 
through a ventilation channel ("vent") and possibly other leaks, form a loop 
which may have a gain higher than 0 dB, at least for some frequencies. With more 
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than 0 dB loop gain, the system may become unstable and produce feedback 
oscillations. 



The present invention is designed to detect a feedback oscillation in the input 
5 signal (55), and set a flag (15) which indicates 'oscillation' or 'no oscillation'. 



Some assumptions about the feedback oscillations in hearing aids are included in 
the design of the detector. The signal produced as a feedback oscillation typically 
consists of a single frequency, namely the frequency at which the loop gain is 
10 highest, taking into account both the linear and non-linear components of the 
hearing aid. The level of the feedback oscillation is relatively stable, after a 
certain settling time. The feedback oscillation often dominates the signal in the 
feedback loop, since its level is often determined by the hearing aid compressors. 



15 The feedback detection process is complicated by the presence of other signals in 
the feedback loop. Many environmental signals, including music, may contain 
segments of periodic nature which may resemble a feedback oscillation. 
However, in the frequency range where oscillations may occur, relatively few 
environmental signals consist of a single frequency only, at least when 

20 considered over a period of a few hundred milliseconds or more. 



The feedback oscillation detector in the present invention is based on measuring 
the overall 'bandwidth' of the signal in the feedback loop consisting of 
components (1), (2), (3), (4), (5) and (50). In the preferred embodiment, the signal 
25 (55) is used as input to the detector, but with slight modifications the detector 

may obtain its input anywhere in the loop. When the bandwidth of the signal (55) 
has been small for a certain minimum period of time, the detector will flag a 
'feedback oscillation' condition. 



30 FIG. 3 describes the detector (49). The input signal (55) is highpass filtered by an 
8-tap FIR filter (36). The filter helps prevent false feedback oscillation detection 
for low frequency input signals since it suppresses the fundamental frequencies 
for a wide range of signals. The 3 dB roll-off frequency for the filter should be 
higher than the lowest expected feedback oscillation frequency. The 8-tap FIR 
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filter is just one example of a usable filter, but many other types may be used. The 
highpass filtered signal (37) is fed to a modeling device (38), which attempts to 
model the spectrum of the signal (37), using a second-order auto-regressive model 
as shown in E4: 

y(n) = x{n) • K - a\y{n - 1) - aiy{n - 2) 
5 (E4) 

where x(n) represents the excitation signal, which drives the model input, while 
y(n) is the output from the model. 

10 The signal model E4 represents a second-order IIR filter with a single complex- 
conjugated pole-pair. Based on the model coefficients a, and a 2 (39), the filters 
center frequency and bandwidth may be computed. This computation is 
performed by the unit (40), which produces a bandwidth (41) and a center 
frequency (48). These two values are compared by (47) to preset thresholds (43) 

15 and (46). The comparator sets flag (44) TRUE if the bandwidth (41) is lower than 
the preset threshold (43) AND the center frequency (48) is higher than the 
acceptable minimum feedback oscillation frequency (46). Otherwise the flag (44) 
is set FALSE. 

20 All components (38), (40), (47) and (45) are working on a frame based schedule. A 
frame length of 40 ms may be used, but other values of the length would also 
work. For each frame, a new value of the flag (44) is computed. Since many 
environmental input signals contain short segments of narrow bandwidth, the 
flag (44) may occasionally be set TRUE while no feedback oscillations are 

25 present. To avoid this, the flag (44) is fed to a stability estimator (45). In here, the 
flag (44) is placed in a delay line which, at any point in time, holds the values of 
the flag from the last N se frames. N se may be selected as 10, but other values would 
also work. The stability estimator (45) sets the detector flag (15) TRUE when and 
only when at least N nin out of the N se past values of the flag (44) were TRUE. For 

30 example, N ln . n may be set to 4. 

The coefficients a, and a 2 in E4 are computed from the autocorrelation 
coefficients R(0), R(l) and R(2), by solving the equations: 
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R(0)-a\ + R(l)-ai = -R(l) 

(E5a) 

R(\)a\ + R(Q)a2 = -R(2) 

(E5b) 

The autocorrelation coefficients may be computed using the following equations: 



(E6a) 



(E6b) 



Hi 



*(1) = ~2>(") *("+1) , n = \..N f -\ 

Nf 



*(2) = ~y>(«)-x(«+2) , «=i-^/-2 

N/ 



(E6c) 



where N r corresponds to the frame length, and x(i) is the i'th sample of signal (37) 
10 from the current frame. 

The 3-dB bandwidth of the filter represented by the auto-regressive model E4 may 
be computed as 

Bandwith = 2(1- Vol) 

(E7) 

15 

and the center frequency may be computed as 



r -1, — "I x 
JCenler = COS ( j=) 



(E8) 



In both equations (E7) and (E8) the result is given in radians. Simple calculations, 
20 in which the system sample rate is included, may be used to convert the values of 



Bandwidth and the f^, into Hz. 
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In the previous description the hearing aid and the methods have been described 
in a simplified manner. Necessary elements like a power source, e.g. a battery, 
and related wiring, the signal processing capabilities of the hearing aid amplifier 
and the interconnecting wiring of the components, as well as the housing, which 
5 is always present have been omitted from the general definition of the hearing aid 
according to the invention. It should be appreciated that these elements of course 
are present in a hearing aid actually manufactured. 
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1. A method for canceling feedback in an acoustic system comprising a 
microphone, a signal path, a speaker, means for detecting presence of 
feedback between the speaker and the microphone and filter means for 
compensating at least partly a possible feedback signal, the method 
comprising: 

- providing a LMS algorithm for generating filter coefficients; 

- where the LMS algorithm operates with a predetermined essentially level 
independent adaptation speed when feedback is not present, this 
representing a first mode, 

where the LMS algorithm operates a level dependent adaptation speed 

when feedback is present, this representing a second mode; 

where the means for detecting the presence of feedback is used to control 

the adaptation mode selection of the LMS algorithm and 

where the feedback detection means comprises bandwidth detection 

means for determining the presence of a feedback signal. 

2. A method according to claim 1, where the update rate for the LMS algorithm 
is determined by the long-term average denominator in the LMS update 
algorithm in the second mode. 

3. A method according to claim 1 or 2, comprising using a highpass filter to 
prevent low-frequency signals from entering the LMS algorithm; where an 
additional feedback cancellation filter and a noise generator is used for 
providing low-frequency input for the LMS algorithm. 

4. A method according to claim 1, where the stability of the signal determined as 
a feedback signal is analyzed. 

5. A method according to claim 4, where the feedback analyzing comprises 
holding flag values from a number of succeeding time frames and comparing 
of these. 
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6. A hearing aid comprising : 
a microphone; 

- a signal path; 
a amplifier; 

5 a speaker; 

- means for detecting feedback between the speaker and the microphone; 

- filter means for compensating at least partly a possible feedback signal; 

- memory means including a LMS algorithm; 

- means for shifting the adaptation mode of the LMS algorithm when 
10 feedback is detected, said means being controlled by the means for 

detecting feedback and 

- means for updating the LMS algorithm by the long term denominator in 
the LMS algorithm; 

where the feedback detection means comprises bandwidth detection 
15 means for determining the presence of a feedback signal. 



7. A hearing aid according to claim 6, comprising stability detecting means for 
the feedback signal. 
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Fig. 1 
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Fig. 2 
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